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ABSTRACT 


Since a planet's magnetic field can sometimes provide a spacecraft with 
some protection against cosmic ray and solar flare particles, it is important 
to be able to quantify this protection. This is done by calculating cutoff 
rigidities. The conventional method of calculating cutoff rigidities in a 
nondipole magnetic field is to trace particle trajectories. This publication 
introduces an alternate method, which is to treat the problem as a boundary 
value problem. In this approach, trajectory tracing is only needed to supply 
boundary conditions. In some special cases, trajectory tracing is not needed 
at all because the problem can be solved analytically. A differential 
equation governing cutoff rigidities is derived for static magnetic fields. 

The presence of solid objects, which can block a trajectory, and other force 
fields are not included. A few qualitative comments, on existence and 
uniqueness of solutions, have been made which may be helpful when deciding how 
the boundary conditions should be set up. The differential equation is first 
expressed in terms of arbitrary coordinates (using the del operator) and then 
it is expressed in terms of a specific set of coordinates (two sets of 
spherical coordinates) . This publication also includes topics on axially 
symmetric fields. If the magnetic field is symmetric about an axis, and if a 
certain kind of cutoff (which this publication calls the "generalized Stormer 
cutoff") is being investigated, the equations simplify considerably. It is 
shown that the vertical cutoff, which is the generalized Stormer cutoff 
evaluated in a direction perpendicular to magnetic east, is constant on a 
magnetic field line. A method for obtaining analytic solutions for a large 
class of problems with axial symmetry is derived and two specific examples are 
included. A procedure for finding the field that provides the best protection 
subject to a certain kind of constraint is also derived and an example is 
given. Each application of this optimization analysis will produce an 
inequality which can be used as an upper bound estimate of the protection that 
can be provided by a given axially symmetric field. 
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Introduction 


Cutoff rigidities are the quantitative measure of the protection provided 
to a region of space by a magnetic field against energetic charged particles. 
This protection is important to a spacecraft in a planet's magnetosphere and 
it also has certain laboratory applications. 

The conventional method of calculating cutoff rigidities in a nondipole 
magnetic field is to trace particle trajectories by computer. There are some 
difficulties (described in Ref. 1) with these methods. This publication 
introduces an alternate method, which is to treat the problem as a boundary 
value problem. A differential equation, which will be called the field 
equation, governing cutoff rigidities for a particle in a static magnetic 
field will be derived. The presence of solid objects, which can block a 
trajectory, and other force fields are not included. If the magnetic field 
has no symmetries, it is not immediately obvious that the field equation 
derived here is more convenient for numerical work than trajectory tracing 
(especially since boundary conditions have to be supplied and they will 
probably come from trajectory tracing) but it does at least provide 
investigators with a choice of methods and they can use their own preferences 
to make a selection. The boundary value problem approach may also offer new 
possibilities for investigating analytic properties of cutoff rigidities. 

In a given magnetic field, it typically happens that there are a number 
of different kinds of cutoff rigidities. For a given location and a given 
particle arrival direction, it is typical to be able to group rigidities into 
discrete bands such that for any rigidity inside one of these bands the 
particle cannot escape to infinity and for any rigidity between the bands the 
particle can escape to infinity (see Ref. 1). In this publication, a "cutoff 
rigidity” refers to the least upper bound or greatest lower bound of any of 
these bands. The field equation will be applicable to any cutoff rigidity 
providing it is well defined everywhere and a differentiable function of its 
arguments. Being well defined means that we cannot investigate a cutoff 
rigidity that "disappears," when position and/or direction are varied, due to 
several distinct bands merging together into one band. In most applications, 
the cutoff of interest will be either the maximum or the minimum of the 
various cutoffs and we expect these cutoffs to be well defined. 
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For a large class of magnetic fields having axial symmetry, it is 
possible to analytically solve the equation for the "vertical cutoff" (the 
cutoff in a direction perpendicular to magnetic east). The solution is first 
derived in the general case and then it is applied to two specific examples. 
One example is the dipole field. The solution is already well known but this 
example is useful for demonstrating the method. The second example is a 
dipole field superimposed on a uniform field in such a way as to create a 
spherical magnetopause. Axially symmetric, isotropically protective fields 
will be defined and these fields will be used when solving for the optimum 
field (the field that produces the greatest magnetic protection) subject to a 
certain kind of constraint. A by-product of the optimization analysis is an 
inequality which can be used as an upper bound estimate of the protection that 
can be provided by a given axially symmetric field. An example will be given. 

This publication assumes that the reader is familiar with the basic 
physical concepts that are discussed in the first few chapters of Rossi and 
Olbert (Ref. 2). Other than that, this publication is a self-contained 
mathematical analysis. 

2. The General Field Equation 


■> 

Let P denote the rigidity of a particle, which is the momentum divided by 
the charge, and let v be the particle velocity. The force equation in a pure 
static magnetic field is 


d_ 

dt 


-* 

(mv) 


qv x B 


where B is the magnetic field in the MKS system of units and q is the charge 
of the particle. The relativistic mass, designated by m, is a constant of 
motion since the speed v of the particle is a constant of motion. It is easy 
to show that the above equation can be expressed in terms of rigidity according 
to 


dP 

dS 


P x B 


( 1 ) 


2 



A ^ 

where P is the unit vector in the direction of P and S is arclength with the 
sense of increasing arclength taken to be in the direction of P. 

Consider a point of observation X and an arbitrary direction represented 
by the unit vector P. Let F(X,P) represent the cutoff rigidity at the point X 
and direction P. Note that there will typically be several cutoff rigidities, 
associated with upper and lower bounds of various rigidity bands. The one 
being investigated is represented by F. Now suppose a particle starts at the 
point X and moves in the direction P and travels with a rigidity equal to the 
cutoff F. The particle is in a state of being on the verge of being able to 
escape to infinity. An infinitesimal change in the rigidity puts the particle 
in a bound state or in an escape trajectory depending on the direction 
of the change. Obviously, as the particle moves to a new point on the 
trajectory, it is still in a state of being on the verge of being able to 
escape to infinity. Therefore, the rigidity of the particle at this new point 
is equal to the cutoff rigidity evaluated at that point and in the direction 
of the particle motion, if the particle has a positive charge (more generally, 
the P in the argument of F refers to the direction of the particle's rigidity 
rather than the direction of motion). Since the particle’s rigidity is 
constant on the trajectory we have 

F(X,P) = constant (2) 

when evaluated on the trajectory of a particle that is initially at some 

* 

point Xq and with initial rigidity in some arbitrary direction P^ with 
magnitude equal to F(X^,P^). The above equation gives 

^ F(X(S) ,P(S) ) = 0 (3) 

where S represents arclength along the trajectory and this trajectory is 

* 

represented by X = X(S) . P(S) is tangent to the trajectory in the direction of 
the particle rigidity. The sense of increasing arclength will be taken to 
be in the direction of P so (1) applies and we also have 
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(4) 


£ = f • 


If we let and be the del operators in position and rigidity space 
and apply the chain rule to (3) and use (1) and (4) we get 


P • V F + (P x B) • V F| =0 

X P lp=F(i,P) 


Equation (5) is an unconventional kind of differential equation because the 

point of evaluation of V^F has P (the magnitude of the particle rigidity) 

equal to the unknown cutoff rigidity. This undesirable characteristic can be 

eliminated the following way. Let P, 0 , <J> denote spherical coordinates in 

P P 

rigidity space. Since F does not depend on the magnitude P, we can write 


■* _ i .3F ; 1 9F ? . 

V p F = i ( 0 p + 4>p) 

F p ae v sin e a<t> v 
p p p 


where 0^ and denote the obvious unit vectors. The above equation applies 

for arbitrary P (note that while F does not depend on P, the above equation 
shows that V^F does). In particular, the equation applies when the derivative 
is evaluated at P = F which gives 


Vlp-F = 1 <— 

V P=F F 30. 


sin 0 V 


The quantity in parentheses was to be evaluated at P = F but since F does not 
depend on P it can be evaluated at arbitrary P. Combining the two equations 
gives 


P 

= - V F 
F p 
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so that (5) becomes 


A -> * -) P -> 

P • V F + (P x B) • - V F = 0 (6) 

x F p 

A *> 

where P is now arbitrary. P and X are also arbitrary because the initial 
conditions of the particle trajectory can be chosen arbitrarily. So in 

4 4 A 

(6), each component of X and of P = PP can be treated as independent variables. 

Equation (6) is one way to write the general field equation. Note that since F 

does not depend on P, V F is perpendicular to P, i.e., 

P 

P • V p F = 0 . (7) 

Equations (6) and (7) and a little vector manipulating produce an alternate 
form 


P 


-» -* p -» 

[V F + B x „ V F] = 0. 
x F p 


( 8 ) 


3. Boundary Conditions, Uniqueness, and Existence of Solutions 

The field equation (equation (6) or (8)) must be supplemented by boundary 
conditions before a solution can be obtained. Existence and, to a lesser 
extent, uniqueness of solutions are sometimes taken for granted in physical 
problems . But existence and uniqueness theorems are still important when 
dealing with boundary value problems because they indicate what kind of 
information and how much information must be supplied in the form of boundary 
conditions before the equation can be solved. This section will not go into 
any elaborate mathematical detail on the subject but will, instead, only make 
a few qualitative comments. 

If equation (1), the force equation, is taken as a background hypothesis, 
the steps that produced (6) from (3) are reversable, so the two equations are 
equivalent. Therefore, existence and uniqueness of solutions to (6) can be 
investigated by investigating those properties of (3). It will be assumed 
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that the boundary conditions have the following basic structure. Suppose we 
are seeking a solution that is to apply to some region R of the 
three-dimensional position space. Let A be the boundary surface of R. It 
might be that A is a closed surface, so that R is a finite region, but not 
necessarily. At any given point on A, we assume that the cutoff has been 
specified in all directions. This is not the only structure that boundary 
conditions can have and from the point of view of existence and uniqueness 
theorems it may not be the best structure. But it is conceptually simple and 
because of this conceptual simplicity, this kind of boundary condition might 
show up frequently in practical applications (note that if the boundary 
conditions are obtained by trajectory tracing, the investigator can make the 
boundary conditions have any structure that he wants). 

It is not difficult to see that existence of solutions will place some 
rather severe restrictions on the boundary conditions that can be used. One 
restriction is the following. Consider a point X on A and an arbitrary 
direction P and let F(X,P) be the cutoff that was supplied as a boundary 
value. Let a particle have initial position X and initial rigidity in the 
direction P and with magnitude equal to the cutoff F(X,P). According to (3), 
the cutoff rigidity is constant on the particle trajectory when evaluated in 
the direction of the particle rigidity (which is tangent to the trajectory). 

If there is another point where the trajectory intersects the surface A and if 

-> ..... 
the point is represented by X* and the direction of the particle rigidity is 

represented by P' , the boundary value at X*,P* must be the same as it was at 
-» - 

X,P, i.e., 


F(X* ,P‘) = F(X,P) . 

The above condition is a special case of a slightly stronger restriction on 

the boundary values, which is as follows. Consider one point X^ on A and a 

direction and another point on A (X^ ^ X^) with another direction P^,. 

Let a particle have initial position X^ and direction P^ and rigidity = 

F(i^,P^) and let another particle have initial position and direction P^ 

with rigidity P„ = F(X ,P ). The condition that must be satisfied is that if 
2 2 2 

P^ ^ P^, then their values must compare in such a way that the trajectories 
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of the two particles do not have a point of tangency. (Note that if P = P , 

12 

a point of tangency is possible because and X^ could be different points on 

the same trajectory so that the laws of mechanics are not violated.) Stated 

another way, if the trajectories have a point of tangency, then P = P . This 

4 12 

can be stated in still another way. Let X^ be an arbitrary point in R and let 
P Q be an arbitrary direction. Consider a set of particle trajectories 
associated with different values of rigidity P. Existence of a solution 
requires that there is not more than one value of P which will produce a 
trajectory such that the boundary value of the cutoff at the intersection of 
the trajectory with A (and in the direction of the trajectory) is equal to P. 
This is a strong requirement on the boundary values and it is not obvious how 
we can test the boundary values to see if the requirement is satisfied. One 
way around this difficulty is to work only with boundary values that represent 
a realistic problem. Boundary values that represent real trajectories in a 
given field will automatically satisfy this requirement. 

Assuming that the conditions required for existence of solutions are 
satisfied, the next question is uniqueness. No matter how small the surface A 
is chosen, uniqueness is guaranteed for at least some arguments of F. F is 

, ,4 a 4 

unique at any point X and direction P such that there exists a point X* on A 
and a direction P' such that a trajectory with initial position X* and 
direction P* and rigidity P = F(X*,P') will reach the point X and direction P. 
There are obviously some choices of X and P that will satisfy this condition. 

4 A 

But for some selected X and P, the requirement that this point and direction 
can be reached by such a trajectory puts some restrictions on where the 

4 A 

surface A can be located. At a point X in the region R and for a direction P, 

4 A 

let F(X,P) be the correct value for the cutoff. The boundary values on the 

4 A 

surface A can only tell us what this value F(X,P) is if the trajectory, 

, 4 A 4 A 

characterized by X,P and P = F(X,P), intersects A. The trajectory is 
undefined at time equal to infinity (a particle with rigidity exactly equal to 
the cutoff is as much bound as it is free) but over any finite distance of 
travel the trajectory is a continuous function of the rigidity. This means 
that over any fixed but finite distance of travel, the trajectory is 
approximately that of a particle that is bound but almost able to escape. 
Therefore, over any finite distance, the particle trajectory is confined to 
the regions of space where the magnetic field is strong enough to keep it 
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bound. In order for the particle to reach the surface A, the surface will 
have to extend into these regions of space. We should anticipate problems 
(i.e., not having enough information to solve the equation) if the boundary 
surface is taken to be an infinitely large sphere or if the magnetic field is 
localized and the surface is taken to be entirely outside of the magnetic 
field. 

In an actual application, the decision on what boundary conditions to 
construct might be made the following way. We might let uniqueness conditions 
decide what boundary surface to use. We would look for a surface (for 
example, an infinite plane) that we expect, through physical intuition, to 
have the property that it will be intersected at least once by the trajectory 
of any particle that is trapped but almost free. If our physical intuition is 
wrong, we might find that there is a limited region in the five-dimensional 
X,P space where the equation can not be solved. Hopefully, this will usually 
not be a serious problem. The boundary values must comply with the existence 
conditions which puts strong constraints on them. One obvious way to comply 
with the existence conditions is to carefully construct the boundary values so 
that they realistically represent the magnetic field. Arbitrarily chosen 
boundary values that are intended to represent a "hypothetical problem" are 
likely to lead to problems. 

The solution to (6) depends strongly on the boundary conditions in the 
sense that it is difficult to extract much information from the equation 
without being specific about the boundary conditions and actually solving the 
equation. We can intuitively see why this is so. The equation only "sees” 
the magnetic field in some limited region of space R. The magnetic field 
could be anything outside of R and the only way the equation "knows" what the 
magnetic field looks like outside of R is through the boundary conditions. It 
would be nice if the region R could be taken to be all space (so that the 
equation will "see" the magnetic field everywhere) and the boundary conditions 
taken to be that F = 0 at infinity. We could probably get more information 
out of the equation (without actually solving it) if we could do this. But as 
already mentioned, there may be some mathematical difficulties with this. 
However, in some special cases (axially symmetric magnetic field) it is 
possible to supplement (6) with a little more physical information and the 
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resulting equation can make some useful predictions without the need of being 
specific about the boundary conditions. Furthermore, less information in the 
form of boundary conditions is needed to solve the equation. This will be 
discussed in a later section. 

4 . Mixed Coordinates 


The position and rigidity gradients in (6) treat X and P as independent 
quantities. However, even though P is not treated as a function of X, the 
coordinates that are selected to measure P may depend on X as well as on P 
(examples are azimuth and zenith angles). Coordinates that are intended to 
refer to P but that actually depend on X as well as on P will be called mixed 
coordinates. If we select a set of coordinates with which to express the 
field equation, we would want the partial derivative with respect to a given 
coordinate to mean that the remaining coordinates are held fixed while the 
given coordinate is varied. does not have that interpretation if mixed 
coordinates are used. This operator varies the position coordinate while 
holding the vector P fixed and this will cause the rigidity coordinates to 
vary. It will be useful to rewrite (6) in terms of different kinds of 
derivatives which have the desired interpretation. 


4 4 4 4 

Let u^X.P) and u 2 (X,P) be two coordinates that are intended to represent 
P. Let us change the function F in (6) so that it refers to the new arguments 

4 4 4 

u and u 2> i.e., so that the cutoff rigidity can be expressed as FCX.u^X.P), 
u 2 (X,P)). Note 



F(X,u 1 (X,P), u 2 (X,P)) 



+ F 


u. 



+ F 


u„ 



where subscripts to the F indicate derivatives with respect to the indicated 
arguments with the other arguments held fixed. Subscripts to F should not be 
confused with components of a vector which are also represented by subscripts. 
In vector form the equation becomes 


-» -» •» -» 

V F = 3 F + F V u, + F Fu_ 
x x u^ x 1 u 2 x 2 
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where 3 denotes a "partial gradient." This is a gradient that holds u and 
constant while varying X. Equation (6) becomes 

P • 3 F + [P • V uj F + [P • V u„] F +(PxB)»~ V F = 0 . (9) 

x xlUj^ x2u 2 F p 

Once the coordinates u^ and u 2 have been selected, expressions for V u^, V X U 2 

and V F can be explicitly calculated and substituted into (9). The resulting 
P 

equation will be in a form that is useful for the coordinates selected. 

Example: Two Sets of Spherical Coordinates 

Consider a coordinate system with rectangular coordinates X,Y,Z and 

AAA 

associated unit vectors e 1 > e 2 ,e 3 ' Let r » e * < t > the usual spherical coordinates 
in that system with unit vectors r,0,<p. Now consider a second coordinate 

AAA 

system X'.Y'.Z' having rectangular unit vectors e’,e’,e' and let this second 

X 4 J 

A A A * 

system be oriented so that e^ is parallel to <f> and e| is parallel to e^. Let 
P,Y»<* be the spherical coordinates in this second system with unit vectors 

AAA 

P,Y,a. A diagram of the coordinate systems together with a list of 
miscellaneous vector identities can be found in Appendix 1. The coordinates 
r,e,<(> will be used to measure X and P,Y.<* will be used to measure P so that 
u^ and u 2 in (9) are given by u^ = a, u 2 = y. 

One of the identities in the appendix is 

a a 

P = tsin y sin a cos $ - cos y sin <(»] e^ 

A A 

+ [sin y sin a sin <t> + cos y cos <|>] e 2 + sin y cos a e 3 (10) 

which gives 


sin y sin a cos 4> - cos y sin <J> = P * e^ 

A A 

sin y sin a sin <(> + cos y cos <j> = P » e 2 

A 

sin y cos a a p • e . 


(ID 

( 12 ) 

(13) 
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In spherical coordinates 


V x ui = V x a = r 


3a 0 3a 
3r + r 30 


+ 


<t! 

r sin 0 


3a 

3<t> 


4 

In these derivatives, P is held fixed while the coordinates r,0,4> are given 
independent variations. Referring to the relative orientation between the 
coordinate systems makes it obvious that if the direction P and the angles 0 
and 4> are held fixed as r is varied, a does not change, which implies 
3a/3r = 0. Therefore, 


■* 0 3a d> 3a 

Vl ' r 30 + r sin 0 3<t> 


Similarly, 


(14) 


3 u = ® ** ± 21 

v x u 2 r 30 r sin 0 d<f> 


(15) 


We can solve for the derivatives in (14) and (15) by differentiating (11), 
(12), and (13) with respect to 0 and <t> and solving the simultaneous equations. 
The equations will not be independent (six equations into four unknowns) so 

A 

we will work only with (12) and (13). Differentiating with P held fixed 
gives 


cos a cos 


cos a cos 


81 

30 


81 

3 <)> 


da 

sin a sin y ~ = 0 
3a 

sin a sm y ~ = 0 

d<J> 


[sin a sin 4> cos 


[sin a sin d> cos 


9 y 3a 

y - cos 4> sin y] T~ + sin y sin 4> cos a ~ 

o<p o<|> 

= cos y sin 4> - sin y sin a cos <t> 

9y 3a 

y - cos <[> sin y] ^ + sm y sin <t> cos a ~ 


= 0 . 


11 


The solution is 


3y _ 3a 
30 " 36 


= 0 


3a 

— = cos a cot y 

dtp 


3<t> 


= sin a 


which gives 


-» 

7 u = 


JL 


x 1 r sin 0 


cos a cot y 


-* 

7 u„ = 


_sL 


x 2 r sin 0 


sin a 


Another identity from Appendix 1 is 


P = sin y cos (a - 0) c + sin y sin (a - 0) 0 + cos y 4> 


(16) 


so that 


2 

p . 3 u = COS a cps__Y 
x 1 r sin y sin 0 


P 



cos y sin a 
r sin 0 


(17) 


Also note 


3 


F = 


x 



0 _ <J> _ 

+ F + , F 

r 0 r sin 0 4> 


which, together with (16), gives 


-» 

3 F = sin y cos (a - 0) 
x 


F + sin y sin (a - 6) p + cos y p 
r r 0 r sin 0 <(> 


(18) 


The derivatives in V^F hold the position coordinates fixed while varying P 

and this means that the primed coordinate system is fixed in space as P is 

varied. Therefore, V can be expressed in terms of the spherical coordinates 

P 

P,Y,a in the usual way, i.e.. 


^ , p 3_ + i 3_ + — a — L. 
p 3P P 3y P sin y d<* 
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But F does not depend on P so we have 


P 3 F = I (yF + .. r g— F ) 
F p F Y sin Y « 


Expressing BasBP+By+Ba gives 
p Y a 


- -> 

PxB = Ba-By 

y ol 


so 


B 


(Px2) • | ^ F “ 1 { - T F - B F } 
F p F sin Y a a Y 


It is more customary to express the magnetic field in terms of the r,0,4> 
components than the P,Y» a components and this conversion can be done by using 
a few more identities in Appendix 1 to get 

^ A ^ * A A 

B = B • Y = B • f cos v cos (a - 6) r + cos Y sin (a - 6) 0 - sin Y $] 

Y 

= B cos y cos (a - 0) + B cos Y sin (a - 0) - B sin y 


and 


B = B • a = B • [sin (0 - a) r + cos (0 - a) 0] 
a 

= B sin (0 - a) + B„ cos (0 - a) 
r 0 
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so that 


A -> P 
(P x B) • : f F 
F p 

1 F a 

= - — [B cos y cos (a - 0) + B Q cos y sin (a - 0) - B, sin y] 

F sin y r <p 


- F [B sin (0 - a) + B cos (0 - a)]} . 
y r 0 


( 19 ) 


Substituting (17), (18), and (19) into (9) finally produces 


... .sin y sin (a - 0) _ . cos y . _ 

[sm y cos (a - 0)] F p + [ r J F 0 + I r sin 0 J ^ 

B B B 

.cos a cos y cot y r , 0 <t», _ 

4- [ ~ — t - Ta + — cot y cos (a - 0) + — cot y sm (a - 0) - zTJ F 

r sin 0 F F Fa 


B B 

.cos y sin a r 0 . _ 

[ — - — + — sm (a - 0) - — cos (a - 0)] F = 0 

r sm OF F y 


( 20 ) 


which is the form of (6) that is appropriate for this particular choice of 
coordinates. 


5. Magnetic Fields With Axial Symmetry 

Other than dipole fields, which are sometimes used to approximate a 
planet’s magnetic field, fields that are symmetric about an axis have limited 
applications in space physics. But investigating such fields may be useful in 
some laboratory applications and it may also be useful in spacecraft 
applications if the idea of using artificial fields to protect small regions 
of space from charged particles is ever taken seriously. Therefore, 
investigating axially symmetric fields seems worthwhile, especially because it 
is possible to go far towards solving the field equation (at least in the case 
of the vertical cutoff) with only a small amount of effort. As already 
mentioned in an earlier topic, when the magnetic field is symmetric about an 
axis, it is possible to supplement the field equation with a little more 
physics so that more information is built into the equation and less 
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information must be supplied in the form of boundary conditions. The analysis 
used here starts with concepts that can be found in Ref. 2. 

We will work with the same coordinates and coordinate systems that were 
described in the example under the discussion of mixed coordinates. These 
coordinates are shown in Appendix 1. The magnetic field is taken to be 
symmetric about the e 3 direction. 

It is shown in Ref. 2 that on an arbitrary particle trajectory there 
exists a constant C that satisfies 


cos y = 


1 

Pr sin 0 


[C 


*<r t 0> 
2ir 1 


( 21 ) 


where P is the particle rigidity (not necessarily equal to the cutoff rigidity) 
and 


$(r,6) = / B • dA (22) 

where the integral is a surface integral on any surface bounded by the circle 
that has the spherical coordinates r and 0. This makes § a function of r 
and 0. C in (21) is constant on a trajectory but it can be different for 
different trajectories. 

Let it be given that a particle was found at a point with coordinates 

r ,0 , y and rigidity P. Solving (21) for C gives 
o o o 

t(r ,6 ) 

C = Pr sin 0 cos y + 

o o 'o 2ir 


Substituting back into (21), we see that a necessary condition for a particle 
to reach a point having coordinates r,0,Y is that the coordinates satisfy 


Pr sin 0 


- [Pr sin 0 cos y + 

4 A A 1 A 


$(r ,0 ) - $(r,0) 


o o 


2n 


-] = cos y 


(23) 
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or 


l .0 ) - *(r,0) 

-1 < r : — - [Pr sin 6 cos y + ] < l 

- Pr sin 0 o o o 2ir — 


which gives 


*(r,0) - t(r ,e ) 

P(r sin 0 + r sin 0 cos y ) > — — 

o o o “ 2ir 


and 


»(r ,0 ) - $(r,0) 

P(r sin 0 - r sin 0 cos y ) > 

o o 'o ~ 2ir 

These inequalities imply: 

If 


$(r,6) > $(r ,0 ) 
o o 


then 


P > 


- 2ir 


$(r,0) - §(r ,0 ) 
o o 

(r sin 0 + r sin 0 cos y ) 
o o 'o 


and 


r sin © + r sin 0 cos y > 0. 
o o o 


If 


$(r,0) > $(r ,0 ) 
o o 
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and 


r sin 0 cos y - r sin 0 > 0 
o o o 


then 


*(r,0) - *(r ,0 ) 

p ^ 9 — 9 

- 2ir(r sin 0 cos y - r sin 0) 
o o o 


If 


$(r,0) < §(r ,0 ) 
o o 


then 


§(r ,0 ) - S(r,0) 

p > 9 — 9 

- 2ir(r sin 0 - r sin 0 cos y ) 
o o o 


and 


r sin 0 - r sin 0 cos y > 0. 
o o o 


If 


*(r,0) < *(r ,0 ) 
o o 


and 


r sin 0 + r sin 0 cos y < 0 
o o o 


then 


*(r ,0 ) - $(r,0) 

p < 9 — 9 

— 2ir|r sin 0 + r sin 0 cos y | 

1 o o o 
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The inequalities that put upper bounds on P indicate the existence of at least 

one band of allowed rigidities where, in this case, "allowed" rigidities 

refer to rigidities such that it is possible to connect the coordinates r Q , 

0 ,y to the point r,0. We will derive a generalization of the Stormer 
o o 

cutoff which is a rigidity that must be exceeded in order for a particle to 
escape to infinity. This kind of cutoff is obtained by working with the 
implications that place a lower bound on P. Those implications are equivalent 
to the statement 

P > f (r ,0 ,y ; r ,0) > 0 (24) 

o o o 

where f is defined by 

§(r,6) - §(r ,6 ) 

o o 

2ir(r sin 0 + r sin 0 cos Y ) 
o o o 

*(r o ,e o ) - $(r ,6) 

2ir(r sin 0 - r sin 0 cos y ) 
o o o 

The condition (24) is a necessary condition for the particle to reach the 

point r,0. (Note that in the arguments that produced (24), it is immaterial 

whether r ,0 ,y are the initial coordinates and r,0 the final coordinates 
o o o 

or vice versa. So f is a cutoff rigidity at r ,0 ,y for particles going to 

o o o 

r,0 or coming from r,0.) A necessary condition for the particle to reach the 
radial distance r is that there exists a 0 that will satisfy (24). But 
there will exist a 0 satisfying (24) if and only if 

P > min f(r ,0 ,y ; r,0) (26) 

0 + ° ° ° 

where the plus sign in the minimization symbol indicates that the minimization 
is done over the values of 0 that result in f being nonnegative, i.e., the 

denominator of the appropriate expression on the right side of (25) must be 
positive. The condition (26) is a necessary condition for the particle to 



if *(r,0) > *(r o ,0 Q ) 

(25) 

if $(r ,0 ) > $(r,0) 
o o 
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reach the radial distance r but there is another, stronger condition. This 

condition is that the particle be able to reach every radial distance r* that 

is between r and r, i.e., we must have 
o 

P > min f(r ,0 ; r f ,0) for every r'e(r ,r) 

- e+ ° ° 0 ° 

which implies 

P > max min f(r ,0 ,y ; r*,0). (27) 

. r , - o o o 

r*e[r ,r] 0+ 
o 

The reader might object to the use of minimums and maximums in (27) instead of 
greatest lower bounds and least upper bounds. This is justified in Appendix 2 
which also discusses some other interesting mathematical properties of the 
cutoff. 

Condition (27) is a necessary condition for the particle to reach the 
radial distance r. A necessary condition for a particle, with initial 
coordinates r ,0 ,y , to escape to infinity is 

Quo 

P > F(r ,0 ) (28) 

o o o 

where we define 

F(r ,0 ) = max min f(r ,0 ,y ; r,0). (29) 

ooo _ * . ooo 

r>r 0+ 
o 

We will call F, defined by (29), the generalized Stormer cutoff. A necessary 
condition for a particle to reach infinity is, according to (28), that its 
rigidity exceed the generalized Stormer cutoff. 

Equations (25) and (29) could be used to solve for the generalized 
Stormer cutoff for a given magnetic field. But these equations are not very 
enlightening in the sense of demonstrating general properties. It is more 
convenient to work with (20) and use (25) and (29) to obtain boundary 
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conditions. However, it is not obvious that the cutoff defined by (29) must 
satisfy (20). The derivation of (20) started with (2). The reasoning that 
produced (2) is intuitively obvious when the cutoff is an upper or lower bound 
to a rigidity band, i.e., when it has the property that an infinitesimal 
change in its value makes the difference between a bound trajectory and an 
escape trajectory. But it wasn’t shown that this is the case for the 
generalized Stormer cutoff. A particle’s rigidity must exceed this cutoff in 
order for the particle to escape but it wasn't shown that if the rigidity does 
exceed this cutoff by a small amount the particle will escape. It is 
therefore not obvious that (2) must apply but it is shown in Appendix 2 that 
(2) does apply and therefore we can use (20). 


Note that the generalization of the Stormer cutoff does not depend on 

a. Assuming that this is the kind of cutoff that is being investigated 

(which will be the case if we are looking for a necessary but not sufficient 

condition for particles initially outside the magnetic field to reach a 

location inside the field), F can be set to zero in (20). Obviously 

a 

F is also zero so the equation, after some rearranging, becomes 
4 > 


6 B 

- . _ sin 6 sin Y _ , r . « . 0 „ , „ 

cos 6 sin Y F r - *■ F 0 - [— sin 0 + — cos 0] F 

r . . „ COS 0 sin Y „ 

= - tan a fsin 0 sin y F + F. 

1 1 r r 0 


+ [=r cos 0 — ~~ sin 0 + co v F } . 
F F r sin 0 y 


But the left side of the above equation does not depend on a and neither does 
the coefficient to the tan a on the right side. So these quantities must 
individually be zero which produces 


F 

r 


tan 0 _ 

F~ 

r 0 


B 

. r tan 0 

l F sin y 


+ 


=— T ] F 

F sin Y Y 


(30) 


F 

r 


+ 


cot 0 
r 



B 

. r cot 0 

1 F s in y 



F sin y 


+ 


COt Y , 

2 1 
r sin 0 



(31) 
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The above equations are the final field equations but it is possible to write 
them in a variety of ways by using one equation to solve for a quantity and 
substituting into the other. One alternate form is 


F + [ COt Y - - — — ] F = 0 
r r F sin y y 


(32) 


rB 

F e + [cot 6 cot T * = 0 


(33) 


Another equation that can be produced from this system is 


(rB + F cot 6 cos y) F + (rB. - F cos y) - = 0 

r r v v v 


(34) 


6. Vertical Cutoff 


The vertical cutoff for a field with axial symmetry is the cutoff in a 
direction that is perpendicular to magnetic east, i.e., y = 90°. Solving 
for the vertical cutoff does not immediately tell us the cutoff for other 
directions but in engineering applications the vertical cutoff is useful to 
the extent that it can give a rough representation of cutoffs for arbitrary 
directions. By assuming that the cutoffs in all directions equal the vertical 
cutoff, we can typically obtain a reasonably accurate approximation of 
directional average fluxes in the case of a magnetic dipole (see Ref. 3) and 
also in the case of using a complex magnetic field which accurately represents 
the Earth’s magnetic field (See Ref. 4. Here vertical cutoff means in the 
direction perpendicular to the surface of the Earth) . We can anticipate that 
this will be the case for many magnetic field patterns and therefore we can 
anticipate that the vertical cutoff will often have useful engineering 
applications. Fortunately, it is easy to solve. 

Setting y equal to 90° in (34) produces 


B e 

Vr + r F e - 0 • 
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This equation is the spherical coordinate representation of the vector equation 


B • VF = 0 

which states that the vertical cutoff is constant along a magnetic field 
line. If the vertical cutoff is known on a surface, it will also be known at 
any point in space that can be connected to that surface by a magnetic field 
line. Therefore, the only thing that requires any work is the construction of 
boundary conditions. 

The boundary conditions are supplied by equations (25) and (29). These 

equations can find the cutoff at any point in space but the computations 

simplify if the equations are used to find the vertical cutoff on the plane 

0 = 90°. At any point in space that can be connected to this plane by a 

o 

magnetic field line, the vertical cutoff can be solved by following the 
magnetic field line until it intersects the plane. When applied to the 
vertical cutoff in the plane, the equations reduce to 

F(r ,6 ,y ) = max min f(r ,90°, 90°; r,0) (35) 

0 ° ° 0 =90°, y =90° r>r 0 ° 

o o o 


where 


f (r ,90°, 90°; r,0) 
o 


|*(r,0) - S(r , 90° ) | 

_ O 

2irr sin 0 


(36) 


Example 1: The Dipole Field 

The Stormer cutoff for a dipole is already well known but the dipole 
provides a simple example to demonstrate the use of the equations in this 
publication. The magnetic field is given by 

B = ~ [2r cos 0+0 sin 0] (37) 

r 


22 


where M is the dipole moment (the dipole moment used here is y Q /4ir times the 
kind of dipole moment that has the dimensions of current times area) . The 

A 

dipole moment vector points in the e 3 direction. The flux is given by 


$(r,6) = JB • dS = 2irr 2 jB r sin 6d6 


4irM .0 . , 

J cos x sin x dx 

r o 


— [1 - cos 26] 
r 


and 


#(r ,90°) 
o 


2irM 

r 

o 


so (36) becomes 


2Mr - Mr + Mr cos 26 

f (r ,90®,90°; r,0) = “ 2 

2r r sin 6 
o 


where we have used the condition r > r to remove the absolute value sign. We 

o 

can see by inspection that the minimum in 0 of F occurs when 0 = 90° and (35) 
becomes 


F <r t 0 n .Y ft ) 
o o o 


M(r - r ) 
o 


e =90°, Y =90 c 
o o 


= max ; 

r>r r 
o o 


Maximizing in r yields 


F(r ,0 ,y ) 
o o o 


M 

4r 2 


e =90°, Y =90* 
o 'o 
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Since only one set of coordinates remain, we will drop the subscripts and write 


F(r,0,Y) 


6=90°, y=90° 



(38) 


which is the vertical cutoff in the equatorial plane. The vertical cutoff is 
constant on a magnetic field line so to find its value at an arbitrary point 
in space, we simply follow the magnetic field line that passes through that 
point to see where it intersects the equatorial plane and use (38). In this 
example it is not difficult to express the solution in equation form. A field 
line satisfies the equation 


dX B 
dS = B ' 


But 


dX dr ~ d© ~ * d* * dr - d© * 
TT = TT r + r ~ 0 + r sin 6„* = Trr + r — 8 
dS dS dS dS dS dS 


so the vector equation produces the two equations 


dr ^r d6 % 

dS ” B ’ r dS B 


which gives 


. B 

dr r 

d6 = r B e 


(39) 


Using (36) gives 


dr 

— = 2r cot 6 
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Separating variables and integrating gives 


r = K sin^ 6 

where K is a constant that identifies the magnetic field line that is being 
investigated. Now consider an arbitrary point r,0. The field line that 
passes through this point intersects the equatorial plane at r^ which satisfies 

r 

o r 

r = — ^ = K = — . 

° sin (90°) sin 0 

2 

The vertical cutoff at r,0 is obtained by evaluating (38) at r = r/sin 0 

o 

and the result is 


F(r,0,9O°) = F(r ,90°, 90°) 
o 


M sin 0 


4r 


(40) 


The well-known solution for the Stormer cutoff for a dipole is (see Ref. 2) 


F(r,0,y) 


H sin^ 0 


2., ,, . 3 ,..1/2,' 

r [1 + (1 - cos y sm 0) ] 


(41) 


which obviously reduces to (40) when y = 90°. In order to predict the more 
general form, equation (41), from the field equation, it would be necessary to 
generalize the boundary values for arbitrary y and then include the y 
dependence by using (32) and (33) or an equivalent system of equations. 

Direct substitution will show that the more general expression given by (41), 
with B given by (37), is a solution to the field equations (32) and (33). 
Incidentally, the dipole field is sufficiently simple that it is fairly 
straightforward to deduce (41) directly from (25) and (29). 
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Example 2: A Spherical Magnetopause 


Let the magnetic field be a superposition of the dipole field in the 

previous example and a uniform field that is antialigned with the dipole 

moment. For some r , the field can be expressed as 

c 


B 


u ~ 2M a 

— [2r cos 6 4- 6 sin 0] - ~ e 
r 6 r^ J 


The physical significance of r can be seen by using the above equation to 

c 

obtain 


B • r = 0 at r = r c . 

4 

In other words, B has no radial component on the sphere of radius r . This 

means that it is permissible to stipulate that the magnetic field is zero for 

r > r , i.e. , the field has a spherical magnetopause (note that the derivation 
c -» 

of (21) made use of B having a zero divergence. This means that a 
magnetopause can’t be chosen arbitrarily because it can't cut across magnetic 
field lines). The magnetic field will be chosen to satisfy 


-» 

B 


M , - * - . .. 2M - 

— [2r cos 0 + 0 sm 0] - “T e, 

J J J 

r 


r < r 


r > r 


The flux is given by 


$(r,0) = 2irM sin 2 0 [- - ^r-] 

r 3 
r 

c 
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for r < r and the flux is zero for r > r . The same methods that have been 
- c c 

used for the dipole can be used for this example and the result is 


r > r 
o c 


F(r ,90°, 90°) 
o 


— [— 
2 l r 
r o 
c 


Ar < r < r 
c — o - c 


(42) 


— { [— 
»r 2 r o 
\ C 


2 

r r 

2 J r* 
r 

c 


— r + — } r < Ar 

,2 r o c 

r’ c 


where 


1/3 


1/3 


A = 


.1/3 


{(1 + /5) + (1 - /5) } 


= (1 + D ) 


2 1/6 {cos (“■ arctan D) - /3 sin arctan D)} 


D = 


r 2 r ^ 

27 t_ 2“ " + 1 

27 r2 r n 


1/2 


Using steps similar to those used in the dipole case, the equation for a 
magnetic field line can be solved and the result is 


r . 2 _ 

~ = K sin 0 

3 3 

r - r 
c 

where K is a constant that identifies the field line. Now let r,0 denote an 

arbitrary point in space. If r > r , the cutoff is zero (the point is outside 

c 

the spherical magnetopause). Otherwise, the cutoff at the point r,0 is given by 

F(r,0,9O°) = F(r ,90°, 90°) 
o 
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where the expression on the right is given by (42) and r Q is the real solution 

(which is unique since r < r ) to the equation 

c 

r 

r o 

, 3 3. .2 3 3 

(r c - r ) sin 0 r c " r o 


7. Isotropically Protective Fields 

In this section we consider fields that are symmetric about the Z axis 
and use the vertical cutoff as the quantitative measure of the protection 
provided by a field. 

Consider a point in space X^. Translate the coordinate system so that 

the point lies in the xy plane and let it have the radial coordinate r . From 

o 

condition (24) we see that the cutoff rigidity for the particle to reach the 

coordinates r ,0 = 90° t y = 90° when it has initial coordinates r,0 is 

oo o 

f (r ,90°, 90°; r,0) . 
o 

This is the protection that the field provides to the coordinates r ,0 = 90° 

o o 

against particles that are initially at r,0. We will say that the field is 

isotropically protective at the point X if for all r > r and 6e(0,90°), the 

o o 

protection is at least as large as the protection against particles that are 
initially at r with 0 = 90°. In other words, the field is isotropically 
protective if 

f(r ,90°, 90°; r,0) > f(r ,90°, 90°; r,90°) 
o o 

for all r > r and all 0c(O,9O°). This condition is, of course, equivalent 
o 

to the statement 

min f(r ,90°, 90°; r,0) = f(r ,90°, 90°; r,90°). (43) 

a o o 
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It was already seen in one of the examples in the previous section that the 
dipole field is isotropically protective, i.e., it satisfies (43), with 
respect to any point in the equatorial plane. 

The practical advantage of isotropically protective fields can be seen 
from equation (35) which implies that the vertical cutoff is bounded according 
to 


F < r .e ,Y ) 

o o o 


< max 
r>r 


f (r , 90° , 90‘ 
o 


r,90°) 


6 =90* 
o 


Y =90' 
'o 


or, using (36), 


F (r .0 ,Y rt ) 
o o o 


< max 
r>r 


|*(r,90°) - *(r ,90°) | 
o 

2irr 


6 =90°, y =90° 
o o 


From (43) we see that the equality is used for isotropically protective fields 

which implies that these fields provide the greatest protection that can be 

obtained for a specified B (r,90°). In other words, if we are given a 

0 

functional form for the 0 component of the field in the xy plane and consider 
the set of all fields that have this 0 component in the xy plane, the greatest 
possible protection is provided by the isotropically protective fields. From 
a computational point of view, these fields have the advantage that the 
vertical cutoff is calculated from the equation 


F(r o’VV 


e =90°, y =90° 
o o 


max 

r>r 

o 


|#(r,90°) - §(r ,90°) | 
o 

2irr 


(44) 


which means that, for this calculation, the field only has to be known in the 
xy plane. Of course something has to be known about the field off of the 
plane in order to conclude that it is an isotropically protective field. 
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There is a systematic procedure for constructing these fields. The first 
step in this procedure is to select a function representing the 0 component 
of the field in the xy plane. This step completely determines the protection 
that will be provided by the field (it also places an upper bound on the 
protection that can be provided by other kinds of fields having the same B 

U 

in the xy plane). The field is not yet uniquely determined. The next step in 

the construction is to select a generating function h(r,0) which will be 

discussed in more detail below. The generating function controls the behavior 

of the field off of the xy plane. In practice, the choice of what generating 

function to use is likely to be determined by the geometry of the region of 

space that is intended to be protected. Since the vertical cutoff is constant 

on a magnetic field line, the protection provided to points in the xy plane at 

the radial distance r is also provided to all points on the surface defined 

by the magnetic field lines that intersect the plane at r = r . The shape of 

o 

this surface can be controlled by making an appropriate selection (by trial 
and error or by other methods) of the generating function. 

We will now investigate where the generating function comes from and how 
it is used to complete the construction of the field. 

The field is to be constructed to satisfy (43) which can be expressed as 

|*(r,6) - $(r o ,90°) | |$(r,90°) - $(r o ,90°)| 

TO ^ n 2irr sin 0 - 2irr 

\7 

The field will be constructed to have mirror symmetry which means that when 
minimizing in 0 we can confine our attention to the interval [0,90°]. A 
sufficient condition to obtain the above equation is that the function 

|*(r,0) - $(r o ,90°) | 

2irr sin 0 

be a decreasing function of 0 on the interval (0,90°). This condition is 
obtained if 


30 


a |*(r,0) - *(r o ,90“)| 
30 2irr sin 0 


for 0e (0 , 90° ) . 


Using 


$(r,0) 


, 2 r e 

2irr J 

o 


B (r,x) sin x dx 
r 


and calculating the derivative gives 


a |*(r,0) - $(r ,90°) | 
d o 

30 2irr sin 0 


l*(r,0) - $(r ,90“) ] 

rB (r,0) - r 2 cos 0 

r n . 2 


2irr sin 0 


if $(r,0) > $(r 


I*(r,0) - $(r ,90°) ] 

-rB (r,0) + cos 0 

r ^ _ . 2 


2irr sin 0 


if $(r ,0) < $(r 


Select a generating function h(r,0) which is arbitrary except for one 
condition which will be stated below. Construct the field to satisfy 

[$(r ,0) - *(r o ,90*)] 

rB (r,0) = h(r,0) cos 0 

2irr sin 0 


The expression for the derivative becomes 


, |*(r,0) - *(r ,90°) | 

o o 

30 2irr sin 0 


l*(r,0) - *(r ,90“)] 

[h(r,0) - 1] 2 cos 0 

2irr sin 0 

if §(r ,0) > *(r 
[$(r,0) - $(r ,90°) ] 

[1 - h(r,0) ] r 2 cos 0 

2irr sin 0 

if $(r,0) < $(r 


(45) 


(46) 


,90“) 


,90“) 


(47) 


,90“) 


,90“) 
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It is evident from the above equation that (45) will be satisfied if 


h(r,0) < 1 for all 0e(O,9O°), r > r . (48) 

“ o 

An isotropically protective field can be constructed from a given function 

B (r,90°) by selecting any function h(r,0) that satisfies (48) and then 
0 

constructing the field from (47). However, (47) is not in a convenient form 
for constructing the field. It is more convenient to work with the auxiliary 
function j which is defined in terms of the generating function by 


tan 0 ~ In j(r,0) = h(r,6) 


(49) 


which can be solved for j giving 


Jh(r,0) cot 0 d0 

J (r,0) = e 


(50) 


The integral in (50) is indefinite because additive constants to the integral 

change j by proportionality constants which have no effect on the final 

results. B can be related to j through (47) and (49) which gives 
r 


-1 


(^ In j (r,0) ) 


sin 0 r 


B r (r,0) 


$(r,0) - $(r ,90°) 
o 

2ir 


Using (46), the above equation can be differentiated to give 

a a - 1 a _1 a 

— l(— In j(r,0)) sin 0] B r (r,0) + [(— In j(r,0)) sin 0] — B r (r,0) 

= B (r,0) sin 0 . 
r 

This equation can be rearranged into 

ii B r <r ' e > - fi lln ‘Ilhii V r - 6) -° ' 
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It is easy to verify that the solution is 


B (r ,0) 
r 


A(r) 3 
sin 0 30 


j(r,0) 


( 51 ) 


where A(r) is an arbitrary function of r. This arbitrariness can be removed by 
the requirement that the field have a zero divergence. A necessary condition 
for a zero divergence is 


„ 2.90° 

2irr J 

o 


B (r,0) sin 0 d0 = $(r,90°) 
r 


(52) 


where 

$(r,90°) = -2ir J r B a (x,90°) x dx. 
O 0 


(53) 


$(r,90°) is regarded as a specified function of r since B (r,90°) is regarded 

0 

as a specified function of r. It is therefore acceptable to have the field 
expressed in terms of $(r,90°). Using (51) and (52) gives 


A(r) 


$(r .90°) 

2irr 2 [j (r,90°) - j(r,0)] 


so (51) becomes 


(04) 


B r (r,0) 


$(r.90°) 

2irr 2 sin 0 Ij(r,90°) - j(r,0)] 


3 _ 

30 


j (r,0) . 


(55) 


Condition (54) is a necessary but not a sufficient condition for the field to 

have a zero divergence. The remaining conditions will be imposed on B (r, 0) 

0 

for 0 * 90°. We require that 


L 3 2 13 

— — (r B (r ,0) ) + : — - — 

2 3r r r sxn 0 30 


(sin 0 B (r,0) ) = 0 
0 
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or, from (55), 


3_ 

ae 


(sin 0 B e (r,6)) = 


1 a_ $(r.90°) 

r 3r zir[j (r,90°) - j(r,0)] 


a_ 

ae 


j (r,6) } . 


Since integration with respect to 0 commutes with differentiation with respect 
to r, the equation can be integrated to give 


sin 6 B (r ,8) 


I , *(r,90 o ) .i (r,6) - .i (r ,0) 
r 3r 1 2ir j(r,90°) - j(r,0) 


(56) 


Equations (50), (55) and (56) completely determine the field components 

in terms of the specified flux function $(r,90°) and the generating function 

h(r,0). If h(r,8) is chosen to satisfy (48) for some r , the field will be 

o 

isotropically protective at the point r = r , 0 = 90°. Altogether, the 

o 

magnetic field is associated with three quantities: The radial coordinate 

(r ) of the point to be protected, the function B (r,90°) (or, equivalently, 
o 0 

the function $(r,90°)), and the generating function h(r,8). In practice it 

may happen that h is constructed so that it explicitly contains r^ to insure 

that (48) is satisfied at the r that has been selected. Whether h does or 

o 

does not explicitly contain the coordinate r , the field will be isotropically 

o 

protective at any point in the plane having a radial coordinate that satisfies 

(48) when substituted for r in (48) . 

o 

Example: Generating function for a dipole 


For an arbitrarily selected r , (47) can be used to solve for the 

o 

generating function from a dipole field associated with the point r = r , 

o 

0 = 90° which is to be protected. The result is 


h(r,0) = 


2 sin 0 

sin^ 0 - — 
r 


It can be seen from inspection that (48) is satisfied, which confirms the 
already known fact that a dipole field is isotropically protective at an 
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arbitrary point in the equatorial plane. This same generating function can 

also be used to construct other fields, associated with other functions 

B (r, 90°), that are isotropically protective at r = r , 0 = 90°. The 
0 o 

auxiliary function j is obtained from (50) and is 


. 2 

r - r sm 0 
j(r,0) f ' 

For a given flux function in the plane, $(r,90°), the components of the field 
are found from (55) and (56) to be given by 


B (r,8) = °?° . s 

irr 


(57) 


sin 0 B Q (r,0) = 


sin 0 3 S(r.90°) . 2. _ . 

- — 5? 2« ■ sin 6 V’ 90 > 


(58) 


Note that a point r was selected to be protected when constructing the 

o 

generating function, but, in this case, the field does not depend on r 

o 

(assuming that B (r,90°) is not chosen to depend on r ) so the field given by 
0 o 

(57) and (58) is isotropically protective everywhere in the equatorial plane. 

If B (r,90°) is chosen to be the function that applies to the dipole field, 

0 

(57) and (58) give the components of the dipole field for 0 ^ 90°. For other 

choices of B (r,90°), other isotropically protective fields will result from 
0 

(57) and (58). A word of caution is in order. The field is required to have 

a zero divergence and this puts a restriction on the choice of functions that 

can be used for B (r,90°). The only functions that can be used are the ones 
0 

that will result in 

$(r,90°) -» 0 as r -» ®. 
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8. Optimum Fields 


This section will investigate the field that provides the greatest 
magnetic protection subject to a certain (admittedly, an unlikely) class of 
constraints. This information could be useful when artificial magnetic fields 
are used for the purpose of providing protection against charged particles. 

By following an appropriate pattern when winding wire around an electromagnet, 
it is possible to create almost any desired field pattern and the question 
arises as to which field pattern is the best. This section answers that 
question for axially symmetric geometries when the constraint involves only 
the normal component of the field in a plane perpendicular to the axis of 
symmetry. Although it is difficult to imagine circumstances where this kind 
of a constraint would be imposed on an experimenter, this section is still 
useful because solving for the protection provided by an optimum field 
provides us with inequalities that can be used to obtain simple upper bound 
estimates of the protection that can be provided by a given field. This is 
demonstrated in the example at the end of this section. 

It is assumed that the set of points where the protection is required is 

a circle of radius r in the xy plane. It is also taken for granted that the 

o 

vertical Stormer cutoff is an adequate description of the protection and that 
will be used as the quantitative measure of the magnetic protection. The 
constraint imposed on the experimenter is assumed to be of the form 

J G(r ,B.(r,90°)) dr = C = constant (59) 

O u 

for some function G which is arbitrary except for a few restrictions which 
will be discussed later. 

The problem of finding the best field subject to the constraint (59) 

comes in two separate parts. The first part is to find the function 

B (r,90°) which maximizes the right side of (44) which is an upper bound on 
© 

the protection that can be provided by any field having the same B (r,90°) 

and is the protection provided by an isotropically protective field. The 

second part is to construct, from this B (r,90°), an isotropically protective 

© 
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field. Such a field is not unique since different generating functions will 
produce different fields. Since the vertical cutoff is constant on a magnetic 
field line, the same protection that is provided to the circle in the xy plane 
is also provided to all points on the surface that is defined by the magnetic 
field lines that intersect the circle. The shape of this surface depends on 
the generating function. If the experimenter has a preference concerning the 
shape of the surface, he might use this as the basis for selecting a 
generating function. 

Methods for constructing isotropically protective fields, when the 

function B (r,90°) is given, have already been provided in the previous 
0 

section. Therefore, in this section we will confine our attention to 

constructing the function B (r, 90°). The objective is to find the 

0 

function that maximizes 

|$(r,90°) - $(r o ,90°) | 

max (60) 

^ 2wr 

r>r 0 

subject to the constraint (59). 

There are a few points concerning constraints that must be discussed. 
Although it is assumed that there is one constraint of the form (59) which can 
be selected arbitrarily (except for a few restrictions), we will actually be 
working with three constraints. The first constraint is somewhat artificial 
and its only purpose is to avoid certain mathematical difficulties. This 
constraint is the requirement that the field be confined to a finite region of 
space, i.e., there exists an A such that the field is zero for r > A. This 
is not a real constraint, in fact it increases the flexibility of the equations 
to follow, because A is arbitrary and we can, if desired, let A go to infinity. 
The primary motivation for imposing this condition is that there are cases 
where the equations to follow will produce mathematically undefined terms if A 
is taken to be infinite during the early steps of the calculations. This 
difficulty can be avoided by taking A to be finite until a final result is 
obtained and then let A go to infinity as the last step. 
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The second constraint is (59) which is supplied by the experimenter. 

Note that imposing some constraint of this form on the field is a necessary 

but not a sufficient condition for the optimum field to exist. To insure that 

an optimum field exists, it will be assumed that G(r,B ) > 0 if B * 0 and 

that G depends on B only through its absolute value, i.e., G(r,-B ) = 

0 0 

G(r ,B 0 ) . 

The last constraint results from the fact that the field has a zero 
divergence. The surface integral of the field on the section of the xy plane 
that lies inside the circle of radius A must be zero. 

For notational brevity, we will write B(r) instead of B (r,90°). It 

0 

should be emphasized that in this section, B(r) represents a component of the 
field rather than the magnitude and it can therefore be negative. The 
constraints are summarized below for easy reference. 


B(r) = 0 for r > A (61) 

J A G(r,B(r))dr = C (62) 

o 

J A rB(r)dr = 0 (63) 

o 

G(r,B) >0 if B * 0 and G(r,-B) = G(r,B) . (64) 

The cutoff is given by 


F(r 0 ) 


max 

re (T q , A ) 


|$(r,90°) - $(r o ,90°) | 
2irr 


where for notational brevity we write F(r ) instead of F(r ,90°, 90°). Using 

o o 

(53) gives 


1 r 

F(r 0 ) = max {“ J pB(p)dp} . 
rc(r 0 ,A) r r Q 


( 65 ) 
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Note that the absolute value sign has been removed in (65). This is because 
we are looking for a field that maximizes the right side of (65) and the 
maximum value that can be obtained is the same whether or not the absolute 
value is retained. A field that produces the maximum in one case can differ, 
at most, by an overall negative sign from a field that produces the maximum in 
the other case and this distinction between the fields is of no consequence. 

The optimum cutoff is the maximum in the function B of the right side of 
(65) which is itself a maximum in r. The maximizing can be done in either 
order and it is convenient to maximize in B first. The optimum cutoff is 
given by 


? . . = max {“ max J [B]} 

optimum rc(r 0 ,A) r r 


( 66 ) 


where the functional J is defined by 

r 

J [BJ = pB(p)dp . 
r r o 


(67) 


We first find the function B that maximizes J subject to the constraints 

r 

(62), (63) and (64). Using standard set theory notation, the set of functions 
that satisfy these constraints can be expressed as 


ae [0,C] {B| *r G( ^* B)d ^ = a > ^ G( * > ’ B)d ' > " C ’ J© pB(p)dp = 0} 


so that 


max J [B] = max U {J [B] | / G(p,B)dp = a, J G(p,B)dp = C, 
_ r r r o 

B ae[0,Cj o 

J A pB ( p ) dp = 0} . 

o 
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The maximum of the union can be obtained by replacing each set in the union 
with the set that contains only the maximum element of the original set. This 
gives 


r 

max J [B] = max U {max{J [B]|J G(p,B)dp = a, 

B F ae[0,C] 

J A G(p,B)dp = C, J A pBdp = 0}} . (68) 

o o 

We now must find the maximum of the set 

G(p,B)dp = a, J A G(p,B)dp = C, J A pBdp = 0} . (69) 

o 

From (67) we see that we are looking for a B that maximizes 

jJJ pBdp (70) 

o 

subject to the constraint 

r 

J r pBdp = - [J pBdp + J* pBdp] (71) 

o 

and the constraints 

J r G(p,B)dp = a (72) 

o 

r o a 

J q G(p,B)dp + J" G(p,B)dp = C - a . (73) 
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To find the maximum of (70), let B, be defined on the interval (r ,r) and be 

* o 

chosen to produce the maximum of (70) subject to the constraint (72). Let 

be defined on (0,r ) U (r,A) and be chosen to produce the maximum of 

o 

r o A 

- [J o pBdp + J r pBdp] (74) 

subject to the constraint (73). There are two possibilities to consider which 
are: 

Case 1: 

assume 


Jp pB l dp - ~ U 0 ° P® 2 dp + pB 2 dpl ' 
o 


(75) 


* 

Then we can find a B, that will satisfy the constraint (73) and produce 


pB^dp = - [J o ° P B*dp + J* pB*dp] 


(76) 


* 

If the strict inequality holds in (75), will not be unique but it is easy 

to see that such a function must exist because it is possible to construct 

such a function out of B^. This can be done by partitioning the domain 

(0,r ) U (r,A) into subintervals and reversing the sign of B on selected 
o 2 

sub intervals so that the various contributions to the integrals in (74) 
subtract with each other instead of add. By condition (64), the constraint 
(73) is still satisfied. Using an appropriate partitioning and sign 
selection, (74) can be given any value between 

r o A 

(J r ° pB 2 dp + J* P B 2 dp] 
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and 


[J o ° pB 2 dp + J* pB 2 dp] 


In particular, it is possible to satisfy (76). Now define 

l B 1 (p) if p e (r Q ,r) 

B*(p) = < * 

I B 2 (p) if p c (0,r Q ) U • 


Note that any element J [B] of the set (69) satisfies (67) subject to the 

r 

constraint (72) and since B^ was chosen to maximize (70) subject to (72) we 
have 


J r [B] < pB^dp = J r [B*] 

o 


i.e., J [B ) is an upper bound for the set (69). Also, since B satisfies (72) 

* r ... * 
and B satisfies (73) and (76), it is evident that B satisfies the 

2 * 
constraints (71), (72) and (73). Therefore J [B ] is an element of the set 

(69) which implies that it is the maximum of the set. To get an explicit 

expression for this maximum value, it is necessary to solve for B^ . The 

maximum of (70) subject to the constraint (72) is a special case of the 

classic prototype problem, treated in the calculus of variations, where the 

integrands may depend on p, B^ and the derivative of B^. The well known 

equation governing the maximizing function is (see Ref. 5) 



d_ 

dp 



) (PB X 


+ X^G) = 0 


where X^ is a constant which is chosen to satisfy the constraint. In our case, 
nothing depends on the derivative of B^ so the equation reduces to 


+ X 


Iff, a 


_a 

dB 


l,r,a 




for p e 


(r ,r) 
o 


(77) 
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The subscripts r and a were added to emphasize the fact that both and X^ 

will in general depend on r and a. This dependence is obtained through (72) 

which must be combined with (77) to produce two simultaneous equations in the 

unknown function B and unknown constant X 

l,r,a l.r.a 

The final conclusion for case 1 is that the maximum value of the set (69) 
is 


J 


r 

r 

o 


pB 


l,r,a 


dp 


where B is obtained from (77) and (72). 

l,r,a 

Case 2: 


Now assume 


pB x dp > - [J o ° pB 2 dp + pB 2 dp] . 


* 

We parallel the steps used in Case 1. Let satisfy (72) 
and 


pB x dp = - [J o ° pB 2 dp + pB 2 dp] 


* 

and define B by 



* 

B^(p) if p e 

B 2 (p) if p e 


<r o’ r) 

(0,r ) U (r, A) 
o 
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Any element J f [B] in the set (69) satisfies 
r o A 

J r [B] = - I / 0 pBAp + / r P Bd <° ] 

subject to the constraint (73) and since B 2 was chosen to maximize (74) 
subject to (73), we have 


J r lB] < - [J o ° pB 2 dp + J* pB 2 dp] = J r [B*] , 


i.e., J [B ] is an upper bound for the set (69). J^[B ] is also an element of 

the set which implies that it is the maximum of the set. To solve for this 

maximum we must solve for B . Note that the conventional derivation of the 

2 

method of Lagrange multipliers applies not only to the case where there is one 

integration interval but also to the case where the region of integration is 

two disjoint intervals. The result is that B^, which will now be represented 

by B to emphasize the dependence on r and a, is given by 

2 ,r ,a 


p + X, 


2 , r , a 8B 


2,r # a 


G (P,B 2> r,a> " 0 


for p e (0,r Q ) U (r.A) 


(78) 


where \ is a constant. The conclusion for Case 2 is that the maximum 

2,r,a 

value of the set (69) is 


- [J ° pB_ dp + pB dp] 

J o H 2 ,r,a H J r ^ 2,r,a r 


where B 0 is obtained from (78) and (73). 
2,r,a 


The two above cases can be combined to give 


max { J r [B] | G(p,B)dp = a, J* G(p,B)dp = C, J* pBdp = 0} 
o 


- min {J r pB l,r,a Ap ' ~ l C pB 2,r,a dp + 4 pB 2,r,a Ap]} 
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Using the above equation together with (66) and (68), it is seen that the 
optimum cutoff compatible with constraints (61), (62), (63) and (64) is given 
by 


optimum re 


, r 

X r o 

max [- max (min { J pB. dp, - If pB dp 
(r ,A) *• ae[0,C] r o 1,r ’ a ° 2 * r ’ a 




(79) 


where B is solved from the system of equations 

l,r,a 


p + X 


l,r,a 3B 


1 , r , a 


G(p,B ) = 0 pc (r ,r) 
l,r,a o 


(80) 


J G(p,B ) dp = a 
r o i,r,a 


(81) 


and B is solved from the system of equations 

2 , r , a 


' + X 2,r,a iT ‘ 0 - c <0 - r o> U (C • A, 

z , r , a 


(82) 


C d e * C 0( 0' B 2.r,a ) ^ * C ' a ' 


(83) 


Equations (80) through (83) will leave B and B undetermined to within 

I* r,a z , r t a 

an overall negative sign. The signs are chosen to make the quantities inside 
the set in (79) positive. This is the sign convention that must be followed 
in order for (80) through (83) to produce maximum values for (70) and (74). 
Maximizing the right side of (79) in a and r can be done using elementary 
methods and is done on a case by case basis. 
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Example: 


In this example the constraint will be that the surface integral of the 
energy density is specified, i.e., 


/ 


A 

o 


pB (p)dp = C 


(84) 


where C is a specified constant. It is difficult to imagine circumstances 
where this constraint would actually be imposed but this example is still 
useful because the results will provide a simple upper bound estimate of the 
protection that can be obtained from a given field. 

Letting 

2 

G(p,B) = pB , 

equations (80) and (82) produce 


B 


l,r,a 


2X 


1 , r , a 


B 


2 , r , a 


2X 


2,r,a 


and (81) and (83) produce 


(85) 


22 2 .2 2 
o r - r , r + A - r 

\ ^ _o 

1 , r , a 8a ’ 2,r,a 8 (C - a) 


The sign of X will be chosen to be negative so that 


J r pB, dp 
r Q r l,r,a 


will be positive and the sign of X will be chosen to be positive so that 

2 ,r,a 


- l C ^.r.a^ + 4 'Vr.a^ 1 
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is positive. Replacing the X*s in (85) gives 


Yr.a - 

o 


1/2 


B 


= _ [ — — ^ 

2 * C * a r 2 + A 2 - r 2 

o 


1/2 


so that 


1/2 


r r _ , .1.2 2., 

'r = l I a (r “ V 1 

o 


1/2 

- .I? 'Yr.a' 1 '’ * ^ Yr.a 4 ' 1 ■ >1 (C - » <r o * * 2 - r2)1 • 


Using (86) and (87), it is easy to show that 


max (min U* P B, dp, - [J ° pB dp + J pB dp]}) 

ac[0,C] r o 1,r ’ a ° Z,r,a r 2,C,a 


2 .2 2 . . 2 2 . 

C(r +A -r)(r - r ) 

. [-* j Y 

2k 


1/2 


where the maximum occurs at 


C , 2 » 2 2 ^ 

a = — (r + A - r ) . 

A 2 ° 


Equation (79) becomes 


l/2\ 


F 


optimum 


max 

rc(r ,A) 
o 


' 2 a 2 2s, 2 2. 

C(r + A -r)(r -r) 

o o_ 

2A 2 


( 86 ) 


(87) 


( 88 ) 
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Maximizing in r gives the final result 




( 89 ) 


(90) 


When a and r are given by (88) and (90), the field is given by 


B(p) 


1/4 



P e 



P e 


(0,r Q ) U 




(91) 


Let us now relax the restriction that the field is confined to a finite 


region 
A -» oo and we have 


The largest value (in A) of F , . occurs, according to (89), when 

optimum 


optimum 

The field that produces this cutoff can be seen, from (91), to be the limiting 
case as the field becomes infinitely "spread out", that is, the field becomes 
infinitely weak but occupies an infinitely large volume of space. This field 
could never be produced but (92) is still useful as an upper bound on the 
protection that can be provided. Consider a field that is axially symmetric 
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but otherwise arbitrary. Assuming that it is a real field (e.g., not the 
ideal dipole), it will have a finite amount of energy stored in it so that 

00 2 

J Q PB dp (93) 

is finite. The protection from this field could not exceed the protection 
from the optimum field subject to the constraint that C is equal to (93) 
evaluated at this specific B, and in fact the protection from this field must 
be less than that from the optimum field (since the optimum field is a 
limiting case that could not exist in reality) so the cutoff for the actual 
field must satisfy the inequality 

1/2 

F < [\ j“ pB 2 dp] . (94) 

Since r^ does not show up on the right side of (94), the upper bound applies 
to all points in the xy plane and also to any point in space that is connected 
to the xy plane by a magnetic field line. 


9 . Summary 

This publication shows that magnetic cutoff rigidities can be regarded 
as a solution to a boundary value problem. This method of solving for the 
cutoff is a possible alternative to trajectory tracing, but in applications, 
trajectory tracing will probably be needed to supply the boundary conditions 
unless the magnetic field has special symmetries. A few qualitative comments, 
on existence and uniqueness of solutions, have been made which may be helpful 
when deciding how the boundary conditions should be set up. The field equation 
can be expressed in the form of (6) or (8). These equations can be applied to 
an arbitrary choice of coordinates providing that the coordinates representing 
the point of observation X and the coordinates representing the direction of 
observation P are separated in the sense that the mapping between the direction 

a -) 

P and the coordinates that represent that direction is independent of X. 

If "mixed” coordinates (which have the property that the coordinates that 
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represent the direction P, depend on X as well as on P) are used, a more 
appropriate form for the field equation is given by (9) which applies to a 
more general class of coordinates. As a special case, (9) was expressed in 
terms of two sets of spherical coordinates. The result was (20). If the 
magnetic field has axial symmetry and the cutoff under investigation is the 
generalized Stormer cutoff, the cutoff is independent of a and satisfies (32) 
and (33). The vertical cutoff, which is the generalized Stormer cutoff 
evaluated in a direction perpendicular to magnetic east, is constant on a 
magnetic field line. This means that the cutoff can be evaluated at an 
arbitrary point in space by following the magnetic field line that passes 
through that point until it intersects a surface where the cutoff has been 
specified as a boundary condition. Analytic expressions for the boundary 
values are given by (35) and (36). Two examples were given to demonstrate the 
methods of calculation. Isotropically protective fields were defined and a 
method for constructing them is given by (50), (55) and (56). These fields 
were used in an optimization analysis which solves for the best magnetic 
protection that can be obtained subject to a constraint of the form (59) . An 
example was given and a by-product of this example is (94) which is an 
inequality that can be used as an upper bound estimate of the protection that 
can be provided by a given field. 
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Appendix 2: Derivation of Equation (2) for the Generalized Stormer Cutoff 


Note that (29) is not a rigorous definition of the Stormer cutoff because 
the existence of maximums and minimums requires that a function and/or its 
domain possess certain properties and it was not shown that this is the case. 

A better definition of the cutoff is 


F <r .rJ = SUP inf f (r .6 ,y ; r,6) 
o o o r>r o o o 

o 6 


(A2.1) 


It will be shown that the cutoff, as defined by (A2.1), satisfies (2) and it 
will also be shown that the definition (A2.1) is equivalent to (29). 


As a reminder, f is given by 


f(r ,0 ,y ; r ,0) 
0 0 0 


*(r,0) - *(r ,0 ) 
o o 

2ir(r sin 0 + r sin 0 cos y ) 
ooo 

*(r ,0 ) - *(r,0) 
o o 

2ir(r sin 0 - r sin 0 cos y ) 
ooo 


if §(r,0) > $(r ,0 ) 

“ o o 

if $(r ,0 ) > *(r,0) 
o o 

(A2.2) 


and the ”+” in the inf symbol in (A2.1) means to include all values of 0 
such that the denominator of the appropriate expression on the right side of 
(A2.2) is positive. 


For arbitrary r ,0 ,y and r > r , let s(r) be the set of values of 0 
ooo o 


satisfying 


0es(r) <=> 


sin 

0 + r 

o 

sin 

0 

o 

cos y 

o 

> 

0 

if 

$(r,0) 

> 

$(r ,0 ) 
o o 

sin 

0 — r 

o 

sin 

0 

o 

cos y 

0 

> 

0 

if 

$(r ,0) 

< 

§(r o* e rt > 

o o 
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so that 


inf f (r ,6 ,y ; r,6) = inf f(r ,6 ,y ; r,6). 
6 + ° ° ° 0es(r) 0 ° 0 

The set s(r) can be expressed as 

s(r) = s^r) U s 2 (r) 

where 

s„(r) = {6 | 6c [0,ir] , *(r,0) > *(r ,6 ), r sin 6 

1 O 0 

+ r sin 6 cos y > 0} 
o o o 


and 


s (r) = {6 | 0e[O,ir], *(r,0) < *(r ,0 ), r sin 0 
2 o o 

- r sin 8 cos y > 0} . 
o o o 

The derivative of f with respect to 8 is easily calculated and i 
below for future reference. 


If (r o’ 0 rt ’^ r ’ e) 

08 000 


r B (r,0) sin 0 
r 


r sin 0 + r sin 0 cos y 
o o o 


[#(r,0) - $(r ,0 )]r cos 0 

O 0 

2ir(r sin 0 + r sin 0 cos y )' 
o o o 


if $(r,e 


Ia (r r,* e rt ’ Y n ;r * 0) 
08 O O O 


-r B (r,0) sin 0 
r 


r sin 0 - r sin 0 cos y 
o o o 


[*(r ,0 ) - $(r,0)]r cos 0 
O o 

2ir(r sin 0 - r sin 0 cos y )' 
0 o o 


if $(r, 


(A2.3) 
(A2 . 4a) 


(A2.4b) 

(A2.4c) 

listed 


) > #<r ,0 ) 
o o 

(A2.5a) 

0) < *<r ,0 n ) 
o o 

(A2.5b) 
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which can also be expressed as 


3f 

ae 


( WV 4 * 0) 



,6 ,y ;r,0) cos 0 - rB (r,0) sin 0] 
o o r 

r sin 0 + r sin 0 cos y 

U ° ° if *(r,0) > #(r o ,0 Q ) 

,0 ,y ;r,0) cos 0 + rB (r,0) sin 0] 
o o r 

r sin 0 - r sin 0 cos y 

° ° ° if #(r,0) < *(r o ,0 o ) 

(A2.6) 


We will now prove several lemmas. The first lemma states that if there 

exists a 6' which is an accumulation point of s(r) (this is a point that can 

be approached arbitrarily closely by elements of s(r)) or is an element of s(r) 

and has the property that §(r,0’) = i(r ,0 ), then the greatest lower 

o o 

bound in (A2.3) is zero. This may at first seem obvious from inspection of 

(A2.2) but it becomes less obvious if we recognize the possibility of 

indeterminate forms. The second lemma describes the geometric structure of 

the set s(r) when r > r and the greatest lower bound in (A2.3) is not 

o 

zero. The third and fourth lemmas will show that the greatest lower bound and 
least upper bound in (A2.1) are a relative minimum and a relative maximum, 
respectively. In other words, F corresponds to a saddle point in f. This 
not only verifies that (29) is equivalent to (A2.1), it also shows (after 
differentiability has been verified) that the extremums are found by setting 
the obvious derivatives equal to zero. This latter result will be important 
in the proof that the cutoff satisfies (2). 


Lemma 1 


Let r ,0 , y be arbitrary coordinates and let r > r and let the magnetic 
o o o o 

field be defined and continuous everywhere. Assume there exists a 0'e[o,ir] 
which is either an accumulation point of s(r) or an element of s(r) and that 
satisfies 

$(r,6’) = $(r ,0 ). (A2.7) 

o o 
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Then 


xnf f(r ,6 ,y ; r,0) = 0. 
n / \ o o o 
6es(r) 


(A2.8) 


Proof 


Inspection of (A2.2) makes (A2.8) obvious if neither r sin 6* + 

r sin 0 cos y nor r sin 0' - r sin 0 cos y are zero. This will be the 
o o o ooo 

case if 0’ is an element of s(r). But since the hypothesis allows for 0' to 
be an accumulation point of s(r) and not necessarily an element of s(r), it is 
possible that one or both expressions are zero. 

First assume that both expresions are zero, i.e., 

r sin 0* + r sin © cos y =0 (A2.9) 

o o o 


and 


r sin 0* - r sin 0 cos y = 0. 

o o o 


(A2.10) 


This implies that r sin 0 cos y =0 and that 0* = 0 or ir. This implies that 

o o o 

$(r,0’) = 0 and from (A2.7) we have $(r ,0 ) = 0. Equation (A2.2) becomes 

o o 

f (r ,0 ,y ;r,0) = 

o’ o’ 'o’ 2irr sin 0 


and (A2.4) gives s(r) = (0,ir). So in this case we have 


inf f(r ,0 ,y ;r,0) = inf l*( r »9Xl 
0cs(r) ooo 0e (0 , ir) 2irr sin e 


Using L’Hopital’s Rule together with 


d$ 2 

“ (r,0) = 2irr B r (r,0) sin 0 


57 


shows that f(r ,0 ,y ;r,0) •* 0 as 0 -* 0 which proves (A2.8) for the special 
o o o 

case that (A2.9) and (A2.10) are both satisfied. 

Now assume that (A2.9) is satisfied but (A2.10) is not satisfied. This 
implies that r sin 0' 4 0. Since the sine is nonnegative on (0,ir) we have 
the strict inequality 

r sin 0* > 0. (A2.ll) 

But (A2.9) and (A2.ll) imply r sin 0 cos y <0 which implies 

o o o 

r sin 0’ - r sin 0 cos y > 0 . (A2.12) 

o o o 

There are two possibilities to consider. The first is that there exists an 
5 > 0 such that either 

0”c(0' ,0* + ?) => $(r,0") < *(r ,0 ) 

o o 


or 


0^(0’ - 5,e f ) => *(r,0") < $(r ,0 ). 

o o 

Assuming that such an £ exists, we can chose it to be sufficiently small such 
that we also have 


r sin 0" - r sin 0 cos y >0. 

o o o 

Then there exists a sequence, {0**}, in s (r) (Equation A2.4) that converges 

n 2 

to 0’. In this sequence we have (using (A2.7) and (A2.12)) 


lim 

n-*» 


$(r,0") - $(r ,0 ) 
n o o 


§(r,0’) - $(r ,0 ) 
o o 


r sin 0’* - r sin 0 cos y r sin 0’ - r sin 0 cos y 
n o o o o o o 


= 0 . 
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This proves (A2.8) when such an 5 exists. Now suppose that no such \ exists. 
Then there exists a sufficiently small £ > 0 such that 

0"e(0 # - 5,0* + V => *(r,0 M ) > $(r ,0 ). (A2.13) 

o o 

But we have the equality when 0” = 0* which implies that $(r ,0 ) is a relative 

o o 

minimum of $(r,0) on the £ neighborhood of 0*. This implies 


8$ 

M <r ’ 6) 



(A2.14) 


Also, from (A2.13), there is a sequence {0"} in s, that converges to 0'. On 

n x 

this sequence we have 


lim f(r , 
o 

n-*» 


0 

o 


;r,0 ) = lim 
n ^ . 

0 -> 0 ' 


*(r,0) - S(r ,0 ) 
o o 

2ir(r sin 0 + r sin 0 cos y ) 
o o o 


To evaluate this limit, note that it is an indeterminate form so use 
L'Hdpital* s Rule to get 


lim f (r ,0 ,y ;r,0") 
o o o n 

n-*» 


lim 

0 -» 0 ' 


31 

30 


(r ,0) 


2ir(r cos 0) 


Note that (A2.9) together with the condition r > r implies that 0* t ir/2. 
Therefore, from (A2.14), the limit is zero. This proves (A2.8) when (A2.9) is 
satisfied. 


Now assume that (A2.10) is satisfied but (A2.9) is not satisfied. Then 

r sin 0' ^ 0 which implies r sin 0' > 0 which implies r sin 0 cos y > 0 and 

o o o 

this implies r sin 0* + r sin 0 cos y >0. Therefore, 0* is an element of 

o o o 

s (r). Obviously, f(r ,0 ,y ;r,0*) = 0 which proves (A2.8). This completes 
1 o o o 

the proof. 
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Lemma 2 


Let r ,0 ,y be arbitrary coordinates and let r > r . Let the magnetic 
o o o o 

field be defined and continuous everywhere. Assume 


inf f (r ,6 ,y ;r,0) > 0 . (A2.15) 

9cs(r) o o o 


Then we have the following implications. 


If: 


r sin 0 cos y =0 (A2.16) 

o o o 


Then: 


s(r) = (0 ,ir) 


(A2.17) 


If: 

r sin 0 cos y >0 and there exists a 0e[O,ir] such that $(r,0) > $(r ,0 ) 
o o o — o o 

(A2.18) 


Then: 


$(r,0) > $(r ,0 ) for all 0e[O,ir] 
o o 


(A2.19) 


and 


s(r) = [0,ir] 


(A2.20) 


If: 


r sin 0 cos y >0 and there exists a 0e[O,ir] 
o o o 


such that §(r,0) < $(r ,0 ) 

~ o o 

(A2.21) 
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Then: 


§(r,0) < §(r o ,e Q ) for all 0elO,ir] (A2.22) 

and 

s(r) = (a,ir-a) (A2.23) 

where a is the solution to r sin a = r sin 0 cos y (A2.24) 

o o o 

If: 

r sin 0 cos y <0 and there exists a 0c[O,w] such that $(r,0) < $(r ,0 ) 
o o o o o 

(A2.25) 

Then: 

$(r,0) < $(r ,0 ) for all 0e[O,ir] (A2.26) 

o o 

and 

s(r) = [0,ir] (A2.27) 

If: 

r sin 0 cos y <0 and there exists a 0c[O,ir] such that #(r,0) > $(r ,0 ) 

O O O “ o o 

(A2.28) 

Then: 

§(r,0) > ^ (r ,0 ) for all 0c[O,ir] (A2.29) 

o o 
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and 


s(r) = (a,ir-a) 


(A2.30) 


where a is the solution to r sin a = -r sin 6 cos y (A2.31) 

o o o 

Proof 


That (A2.17) follows from (A2.16) is obvious from inspection of (A2.4). 

Now assume (A2.18). Since r sin 0 cos y > 0 we have 

o o o 


{0 | 0 e[0,ir], r sin 0 + r sin 0 cos y > 0} = [0,ir] 


and 


{0 I 0 e[0,ir], r sin 0 - r sin 0 cos y > 0} = (a,ir-a) 

o o o 

where "a" is given by (A2.24). s^r) and s^r) (Eq. A2.4) become 

s(r) = {0 | 0 c [0 ,ir] , §(r,0) > $(r ,0 )} (A2.32) 

1 o o 

s (r) = {0 | 0 e(a,ir-a), $(r,0) < $(r ,0 )} (A2.33) 

2 o o 

By hypothesis, there exists a 0’e[O,ir3 such that $(r,0') > fCr^.r^). We will 

show that $(r,0) > §(r ,0 ) for all 0e[O,ir] by contradiction. Assume there 

o o 

exists a 0”e[O,ir] such that #(r,0") < $(r ,0 ). But $(r,0) is continuous so 

o o 

the conditions §(r,0’) > $(r ,0 ) and $(r ,0”) <§(r ,0 ) imply that there 

o o o o 

exists a 0"*elO,ir] such that $(r,0''') = $(r ,0 ). This implies that 0' ' * e 

o o 

s^(r) which contradicts (A2.15) and Lemma 1. This proves (A2.19). But (A2.19) 
implies that s^Cr) = [0,ir] which proves (A2.20). 

Now assume (A2.21). As in the previous case, s^(r) and s^Cr) are given 

by (A2.32) and (A2.33) with ”a” the solution to (A2.24). By hypothesis there 

exists a 0e[O,ir] such that $(r,0) < $(r ,0 ). But this violates the conclusion 

o o 

(A2.19) and therefore the hypothesis (A2.18) must be false, i.e., there is no 
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6c [0 ,«] such that f(r,0) > #(r o ,0 Q ). This proves (A2.22). But (A2.22) implies 
that s^tr) is empty and s 2 (r) = (a,ir-a). This proves (A2.23). 

Now assume (A2.25). The same steps that produced (A2.32) and (A2.33) 
will now produce 

s (r) = {0 | 0 e(a,ir-a), $(r,0) > *(r ,0 )} (A2.34) 

1 o o 

s (r) = {0 I 0 c[0,ir], *(r,0) < *(r ,6 )} (A2.35) 

2 o o 

where "a** is the solution to (A2.31). By hypothesis there exists a 0'e[O,ir] 

such that $(r,0') < $(r ,0 ). We will show that $(r,0) < $(r ,0 ) for all 

o o o o 

0e[O,ir] by contradiction. Assume there exists a 0”e[O,ir] such that §(r,0") 

> §(r ,0 ). But $(r ,0) is continuous so the conditions §(r,0’) < $(r ,0 ) 

~ O o o o 

and $(r,0") > #(r ,0 ) implies that there exists a Q’ ’ ’ with the following 
o o 

properties: 

$(r,0* • •) = #(r ,0 ) (A2.36) 

o o 


and 


Every deleted 5 neighborhood of 0* * ' contains points such that 

*(r,6) < *(r ,6 ) (A2.37) 

o o 

From (A2.37) we have it that O'” is an accumulation point of s 2 (r). Therefore 
(A2.36) and Lemma 1 produce a contradiction. This proves (A2.26). But (A2.26) 
implies that s 2 (r) = [0,ir] which proves (A2.27). 

Now assume (A2.28). But this violates the conclusion (A2.26) so the 
hypothesis (A2.25) must be false, i.e., there is no 0e[O,ir] such that $(r,0) 

< i(r ,0 ). This proves (A2.29). But (A2.29) implies that s (r) = (a,ir-a) and 

O O 1 

s^Cr) is empty. This proves (A2.30). 
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Lemma 3 


Let r ,0 ,y be arbitrary coordinates and let r > r . Let the magnetic 
o o o o 

field be defined and continuous everywhere. Assume 


inf f (r ,6 ,y ;r,0) > 0. (A2.38) 

_ , . o o o 

6cs(r) 


* 

Then there exists a 0 (r)es(r) such that 


$(r,0*(r) ) * *(r ,0 ) 
o o 


(A2.39) 


inf f(r ,0 ,y ;r,0) = f(r ,0 ,y ;r,0 (r)) (A2.40) 

^ x OOO 000 

0cs(r) 


3f 

30 


(r ,0 ,y ;r,0) 
0 0 0 


= 0 


0 = 0*(r) 


(A2.41) 


Proof 


* 

It suffices to show that there exists a © (r)cs(r) that satisfies (A2.40) 
and (A2.41) since (A2.39) follows from (A2.38) and Lemma 1. 


First consider the special case where r sin 0 cos y =0. Then s(r) = 

000 

(0,ir) and 


f(r o’VV r ’ e) = 


|$(r,0) - Hr ,0 ) | 
o o 

2irr sin 0 


From (A2.38) and Lemma 1 we have it that $(r,0) ^ Hr ,0 ) and $(r,ir) ^ 

o o 

§(r ,0 ). Therefore f increases without bound as 0 -* 0 or 0 -* ir. The greatest 
o o 

lower bound must be a relative minimum at an interior point in (0,ir) which 
proves the lemma for this special case. 
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Now assume 


r sin 0 cos y ^ 0. 
o o o 

Lemma 2 states that s(r) is an interval. If it can be shown that the greatest 

lower bound is not obtained by taking the limit as 0 approaches an end point 

of the interval it will follow that the greatest lower bound is a relative 

minimum which occurs at some interior point in the interval. This will imply 

* 

that there exists a 0 (r)cs(r) satisfying (A2.40). It will also imply that 
* 

0 (r) satisfies (A2.41) providing the derivative exists. Inspection of (A2.2) 
and the possible intervals listed in Lemma 2 shows that under each possible set 
of conditions ((A2.18), (A2.21), (A2.25) and (A2.28)), f is differentiable on 
the interior of the interval that applies to those conditions. Therefore, 
(A2.40) and (A2.41) can both be proven by showing that 

If: 


r sin 0 cos y > 0 and $(r,0) > *(r ,0 ) for all 0c[O,ir] (A2.42) 

o o o oo 


Then: 


lim f(r ,0 ,y ;r,0) > inf 
0->O 0 0 0 0ctO,ir] 


f (r ,© ,Y ;r,0) 
o o o 


(A2.43) 


and 


lim f(r ,0 ;r,0) > inf f(r ,0 t y ;r,0) (A2.44) 

* ^ OOO _ r . , ooo 

0-»ir 0c[O,ir] 


If: 


r sin 0 cos y >0 and §(r,0) < $(r ,0 ) for all 0e[O,ir] (A2.45) 

ooo o o 


Then: 


lim f(r ,0 ,y ;r,0) > inf f(r ,0 ;r,0) (A2.46) 

„ OOO - . . ooo 

0-»a 0c(a,ir-a) 
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and 


lim f(r ,6 ,y ;r,0) > inf 
0-»ir-a o o o 6e(a,ir-a) 


f (r ,0 ,y ;r,0) 
o o o 


where 


r sin a = r sin 0 cos y 
o o o 


If: 


r sin 0 cos y <0 and $(r,0) < $(r ,0 ) for all 0e[O,ir] 
o o o o o 


Then: 


lim f(r ,0 ,y ;r,0) > inf 
6-*0 o o o 0e [0,ir] 


f (r ,6 .Y rt ; r . e ) 

o o o 


and 


lim f(r ,0 ,y ;r,0) > inf 
0-»ir 0c [ 0 , ir ] 


f (r .Y^;r,0) 
o o o 


If: 


r sin 0 cos y <0 and $(r,0) > $(r ,0 ) for all 0e[O,ir] 
o o o o o 


Then: 


lim f(r ,0 , y ;r,0) > inf f(r ,0 ,Y ;r,0) 
. o o o - , N o o o 

0->a 0e(a,ir-a) 


and 


lim f(r ,0 ,y ;r,0) > inf 
^ o o o , x 

0-»ir-a 0c(a,ir-a) 


f (r ,0 ,y 
o o o 


(A2.47) 


(A2.48) 


(A2.49) 


(A2.50) 


(A2.51) 


(A2.52) 


(A2.53) 


(A2.54) 
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where 


r sin 0 = -r sin 0 cos y (A2.55) 

o o o 

Inspection of (A2.2) shows that (A2.46) and (A2.47) follow from (A2.45) because 

the limits produce singularities. Similarly, (A2.53) and (A2.54) follow from 

(A2.52) (note that §(r,a) ^ $(r ,0 ) and §(r,ir-a) 4 $(r ,0 ) because a and ir-a 

o o o o 

are accumulation points of s(r)). Now assume either (A2.42) or (A2.49). Since 

r sin 0 cos y 0, then in either case the derivative of f (given by (A2.5)) is 
o o o 

negative in the limit as 0 -» 0 and positive in the limit as 0 ■* ir. Therefore 
f is decreasing when 0 is in a sufficiently small neighborhood of zero and 
increasing when 0 is in a sufficiently small neighborhood of ir. The least 
upper bound cannot occur in either of these limits which proves (A2.43), 

(A2.44), (A2.50) and (A2.51). This completes the proof. 


Lemma 4 


Let r ,0 ,y be arbitrary coordinates. Let the magnetic field be 
ooo 

differentiable in r,0 for all r > 0. Assume that $(r,0) is bounded as r ■» ®. 
Assume also that 


sup inf f(r ,0 ,y ;r,0) > 0. (A2.56) 

r>r 0cs(r) ooo 
o 


* * * 

Then there exists a finite r > r and a 0 es(r ) that satisfy the 

o o o o 

following conditions: 


sup inf f(r ,0 ,y ;r,0) = f(r ,0 ,y ;r* ,0* ) (A2.57) 

^ /v/x OOO ooooo 

r>r Gcs(r) 
o 


fe <r o'W r - e) 



(A2.58) 


(r ,0 ,y ;r,0) 
3r ooo 



(A2.59) 
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Proof 


From (A2.56) we have it that there exists a set R such that 


rcR => r > r and inf 

° 0es(r) 


f (r ,0 , y ;r,0) > 0 
o o o 


( A2 . 60) 


and 


sup inf f(r ,0 ,y ;r,0) = sup inf f(r ,0 ,y ;r,0). 

r>r 0es(r) o o o rc R 6cs(r) 0 ° ° 

o 


From Lemma 3 we have 


sup inf f(r ,0 , y ;r,0) = sup f(r ,0 ,y ;r,0*(r)) 
^ . ooo _ ooo 

r>r 0es(r) rcR 

o 


(A2.61) 


for some 0*(r)es(r) which satisfies (from (A2.41) and (A2.6)) 
f (r o ,0 o> Yo! r » e *( r )) cos ©*(r) = rB r (r,0*(r)) sin 0*(r) 

if $(r,0*(r)) > $(r o ,0 o ) (A2.62) 

f (r o ,0 o ,Y o ;r,0*(r)) cos 0*(r) = -rB r (r,0*(r) ) sin 0*(r) 

if $(r ,0*(r) ) < $(r o ,0 o ). (A2.63) 

Note that one expression or the other applies by virtue of (A2.39). First 
assume that R is bounded above. Then for r sufficiently large we have 


inf 

0es(r) 


f (r ,© ,Y„; r ,e) 

OOO 


f(r «’ e „’Y -.r,e*(r)) = o 

OOO 


so the least upper bounds in (A2.61) must occur at some finite r . Now 

o 

assume that R is not bounded above. By hypothesis $(r,0) is bounded and this 

implies that rB (r,0) -» 0 as r -» «. From (A2.62) and (A2.63) we have it 
** * * 

that f(r ,0 ,y ;r,0 (r)) -» 0 as r -» » unless 0 (r) -» ir/2. But from (A2.2) we 
o o o . 

see that f (r o ,0 Q ,Y o ;r,0 (r)) still goes to zero as r ■> ® even if © (r) does 
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* 

approach ir/2. Therefore, the least upper bound of f(r ,0 ,y ;r,0 (r)) must 

jt o o O 

occur at some finite r in all cases. Note that we must have r > r 
* ° o o 

because if r 0 = r Q > f can t> e made equal to zero by letting 0 = 0 q and this 

implies that f(r ,0 ,y ;r ,0 (r )) = 0 which is clearly not the least upper 
0 o o o o 

bound in (A2.61). This proves that 


sup inf f(r ,0 ,y ;r,0) 
o o o 

r>r 0cs(r) 
o 


f(r 0 ,y ;r* 0*(r*)) 
O O 0 o o 


(A2.64) 


* 

for some finite r > r . 

o o 


This proves (A2.57). Note that 


f(r ,0 ,y ;r* 0*(r*)) = inf f(r ,0 , Y :r* 0) 
oooo o . oooo 

0cs(r*> 

0 


(A2.65) 


Conditions (A2.56), (A2.64) and (A2.65) imply that the hypothesis to Lemma 3 

* 

is satisfied at r = r so (A2.39) gives 

o 

*(r , 0 (r )) * *(r ,0 ) . 
o o o o 

* * 

This means that for r in a sufficiently small neighborhood of r Q , $(r,0 (r)) - 

§(r ,0 ) does not change sign, so f(r ,0 ,y ;r,0 (r)) is a well-behaved function 
o o o o o * 

of r in that neighborhood, i.e,, f is differentiable at r . Since the least 

° * 

upper bound in (A2.61) is a relative maximum at an interior point r^ and f is 
differentiable there, the derivative is zero there, i.e., 


df * 

a^ (r «» e «*Y rt ;r,0 (r)) 
or o o o 


= 0. 


r = r 
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The chain rule gives 


— “(r ,0 ,y ;r,0) 
dr o’ o 'o’ * 


* 

r 

o 

* * 
0 (r ) 


+ ^ ( V 6 o*V r * e) 


36* (r) 
3r 


* * 

0 

o 


The above equation together with (A2.41) prove (A2.58) and (A2.59) where 
* * * 

6 = 0 (r ). This completes the proof, 

o o 

We are now in a position to verify (2). Let a particle have initial 

coordinates r ,0 ,y and let the particle move with a rigidity P equal to the 
o o o 

cutoff F(r ,0 ,Y ) which is assumed to be greater than zero. Lemma 4 states 
o o o 

that the rigidity can be expressed as 


P = f (r ,0 ,y ;r*,0*) 
o o o o o 


Let the particle move a small distance dfi. to the new coordinates = r^ + dr, 

0 =0 + d0 . y = y + dy. To first order, the cutoff rigidity at the new 
1 o 10 

coordinates is 


f<r 1 ,e i , V r*.e* ) . f(r 0 .« 0 .T 0 i<,«*) 


3f 

(r ,0 ,y ;r*,0*) 

df 

dr + — - 

(r ,0 

8r 

0 

o o o o o 

30 

o 

O ( 

df 

(r ,0 ,y 

dY + TT 

(r ,0 

3Y 0 

o o o o o 

3r 

o 

0 

df 

* 

30 

o 

( r .e ,Y n ; r ^ e * ) 

o o o o o 

d0* . 



* \ j * 


(A2.66) 
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To be definite, suppose that $(r£,e£) > $(r o ,0 o ). Analogous steps will reach 
the same conclusion if $(r£,0£) < $(r o ,0 o ) (note that the equality is ruled 
out by (A2.39) together with 0 q = 0*(r£)). Then we use the upper expression 
on the right side of (A2.2) when taking the derivatives and we obtain 


3f_ 

3r 


3f 

30 


M. 


3*(r ,0 ) 

o o 

3r 

o 

* k 

2ir(r sin 0 + r sin 0 cos y ) 

o o o o 'o 


3$(r ,0 ) 
o o 


30 

o 

* * 

2ir(r sin 0 + r sin 0 cos y ) 

o o o o o 


[*(r ,0 ) - #(r ,0 ) ] sin 0 cos y 
o o o o o o 

ic "k 2 

2ir(r sin 0 + r sin 0 cos y ) 
o o o o o 

k k 

[§(r ,0 ) - *(r ,0 )]r cos 0 cos y„ 
o o o o o o o 

k k 2 

2ir(r sin 0 + r sin 0 cos y ) 
o o o o o 


k k 

[*<r ,e ) - *(r ,0 ) J r sin 0 sin y 
o o o o o o o 

* * ^ 
2ir(r sin 0 + r sin 0 cos y ) 

o o o o o 


The last two derivates on the right side of (A2.66) are zero because of 

* * 

Lemma 4. Expressing the terms that contain $( r ,0 ) - #(r ,0 ) in terms of 
* * o o o o 

f(r ,0 ,y ;r ,0 ), which in turn is equal to P, and substituting into (A2.66) 
o o o o o 

gives 


k k k k 

f (r ,0 ,y ;r ,0 ) - f(r ,0 ,Y_;r ,0 ) 
11111 ooooo 


k k 

v sin 0 + r sin 0 cos y 

o o o o o 


1 3$ 

{ - [ 7“ 7~ (r ,0 ) + P sin 0 cos y_J dr 
2ir 3r o o o o 

o 


1 3$ 

- [~ TT - (r ,0 ) + Pr cos 0 cos y 1 <1© + Pr si- n © si n T (A2.67) 

2ir 30 o o o o o o o 'o' 

o 


Note that (21) gives 


* ( r o , 0 o ) 

Pr sin 0 cos y + r = constant. 

o o o 2v 

Taking the differential of the above equation shows that the right side of 
(A2.67) is zero, i.e., the cutoff is constant on the trajectory. 
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